home *** CD-ROM | disk | FTP | other *** search
/ CU Amiga Super CD-ROM 19 / CU Amiga Magazine's Super CD-ROM 19 (1998)(EMAP Images)(GB)[!][issue 1998-02].iso / CUCD / Programming / LEDA / source / src / plane / _Point1.c.ORG < prev    next >
Encoding:
Text File  |  1994-11-16  |  3.2 KB  |  160 lines
  1. /*******************************************************************************
  2. +
  3. +  LEDA  3.1c
  4. +
  5. +
  6. +  _Point1.c
  7. +
  8. +
  9. +  Copyright (c) 1994  by  Max-Planck-Institut fuer Informatik
  10. +  Im Stadtwald, 6600 Saarbruecken, FRG     
  11. +  All rights reserved.
  12. + 
  13. *******************************************************************************/
  14.  
  15.  
  16. #include <math.h>
  17. #include <ctype.h>
  18. #include <LEDA/Point1.h>
  19.  
  20. //------------------------------------------------------------------------------
  21. // Points 
  22. //------------------------------------------------------------------------------
  23.  
  24.  
  25. Point_rep::Point_rep()  { count=1; x = y = 0; w = 1;}
  26.  
  27. Point_rep::Point_rep(Int a, Int b) 
  28. { x = a; 
  29.   y = b;
  30.   w = 1; 
  31.   count = 1; 
  32. }
  33.  
  34.  
  35. Point_rep::Point_rep(Int a, Int b, Int c) 
  36. { x = a; 
  37.   y = b;
  38.   w = c; 
  39.   count = 1; 
  40. }
  41.  
  42. Point::Point()                  { PTR = new Point_rep; }
  43. Point::Point(Int x, Int y){ PTR = new Point_rep(x,y); }
  44. Point::Point(Int x, Int y, Int w){ PTR = new Point_rep(x,y,w); }
  45.  
  46. //Point::Point(vector v)          { PTR = new Point_rep(v[0], v[1]); }
  47.  
  48. /*
  49. double Point::angle(const Point& q, const Point& r) const
  50. {
  51.   double cosfi,fi,norm;
  52.   
  53.   double dx  = q.ptr()->x - ptr()->x; 
  54.   double dy  = q.ptr()->y - ptr()->y; 
  55.  
  56.   double dxs = r.ptr()->x - q.ptr()->x; 
  57.   double dys = r.ptr()->y - q.ptr()->y; 
  58.   
  59.   cosfi=dx*dxs+dy*dys;
  60.   
  61.   norm=(dx*dx+dy*dy)*(dxs*dxs+dys*dys);
  62.  
  63.   cosfi /= sqrt( norm );
  64.  
  65.   if (cosfi >=  1.0 ) return 0;
  66.   if (cosfi <= -1.0 ) return M_PI;
  67.   
  68.   fi=acos(cosfi);
  69.  
  70.   if (dx*dys-dy*dxs>0) return fi;
  71.  
  72.   return -fi;
  73. }
  74.   
  75.  
  76. Point Point::rotate(const Point& origin, double alpha) const
  77. { if (origin == *this) return *this;
  78.   segment s(origin,*this);
  79.   return s.rotate(alpha).end();
  80. }
  81.  
  82. Point Point::rotate(double alpha) const
  83. { return rotate(Point(0,0),alpha); }
  84.  
  85. Point Point::translate(double alpha, double d) const
  86. { double dx = cos(alpha) * d;
  87.   double dy = sin(alpha) * d;
  88.   return Point(ptr()->x+dx,ptr()->y+dy);
  89.  }
  90.  
  91. Point Point::translate(const vector& v) const
  92. { return Point(ptr()->x+v[0],ptr()->y+v[1]);
  93.  }
  94.  
  95. double Point::distance(const Point& p)  const
  96. { return hypot(p.ptr()->x - ptr()->x, p.ptr()->y - ptr()->y); } 
  97.  
  98. double Point::distance() const
  99. { return distance(Point(0,0)); }
  100. */
  101.  
  102.  
  103. int Point::operator==(const Point& p) const 
  104. { return ( (((ptr()->x) * (p.ptr()->w) - (p.ptr()->x) * (ptr()->w)) == 0) && 
  105. (((ptr()->y) * (p.ptr()->w) - (p.ptr()->y) * (ptr()->w)) == 0)); }
  106.  
  107.  
  108. ostream& operator<<(ostream& out, const Point& p)
  109. { out << "(" << p.X() << "," << p.Y() << "," << p.W() << ")";
  110.   return out;
  111.  } 
  112.  
  113. istream& operator>>(istream& in, Point& p) 
  114. { // syntax: {(} x {,} y {,} w {)}   
  115.  
  116.   int x,y,w;
  117.   char c;
  118.  
  119.   do in.get(c); while (in && isspace(c));
  120.  
  121.   if (!in) return in;
  122.  
  123.   if (c != '(') in.putback(c);
  124.  
  125.   in >> x;
  126.  
  127.   do in.get(c); while (isspace(c));
  128.   if (c != ',') in.putback(c);
  129.  
  130.   in >> y; 
  131.  
  132.   do in.get(c); while (isspace(c));
  133.   if (c != ',') in.putback(c);
  134.  
  135.   in >> w; 
  136.  
  137.   do in.get(c); while (c == ' ');
  138.   if (c != ')') in.putback(c);
  139.  
  140.   p = Point(x,y,w ); 
  141.   return in; 
  142.  
  143.  } 
  144.  
  145. /*
  146. int compare(const Point& a, const Point& b)
  147. { Int x = a.X() * b.W();
  148.   Int y = b.X() * a.W();
  149.   int s = sign(a.W()) * sign(b.W());
  150.   if (x == y) 
  151.   { x = a.Y() * b.W();
  152.     y = b.Y() * a.W();
  153.     if (x == y) return 0;
  154.    }
  155.   return (x < y) ? -s : s;
  156. }
  157. */
  158.  
  159.  
  160.